Quantum fluids of self-assembled chains of polar molecules 
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We study polar molecules in a stack of strongly confined pancake traps. When dipole moments 
point perpendicular to the planes of the traps and are sufficiently strong, the system is stable 
against collapse but attractive interaction between molecules in different layers leads to the formation 
of dipolar chains, analogously to the chaining phenomenon in classical rheological electro- and 
magnetofluids. We analyze properties of the resulting quantum liquid of dipolar chains and show that 
only the longest chains undergo Bose-Einstein condensation with a strongly reduced condensation 
temperature. We discuss several experimental methods for studying chains of polar molecules. 



Recent progress in trapping and cooling of chromium 
atoms [l[ and polar molecules 0] opened new directions 
for studying quantum systems with dipole interactions. 
Manifestations of dipolar interactions have been observed 
in the density profiles of the chromium condensate^ and 
predicted in the dispersion of the Bogoliubov mode [|| . 
Exotic many-body states arising from the anisotropy and 
long range character of dipolar interactions have also 
been theoretically proposed @@]. An important aspect 
of dipolar systems is the attractive part of the inter- 
action for dipoles aligned head to tail, which can lead 
to the system collapse @- To circumvent this problem, 
most theoretical analysis focused on two and one dimen- 
sional systems with dipolar moments perpendicular to 
the plane(axis) of the sample, in which case the attrac- 
tive part of the interaction is absent. 

In this paper we consider polar molecules in a three 
dimensional trap which was divided into a stack of two- 
dimensional (2D) pancakes by applying a strong one- 
dimensional optical lattice @(see Fig. [Ha)). When dipo- 
lar moments are polarized in the direction perpendicular 
to the planes, the intra-layer interaction is purely repul- 
sive @. The attractive interaction, however, is present 
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FIG. I: (a) Schematic figure of a stack of pancake traps. 
External electric held polarizes molecules in the direction 
perpendicular to the layers. Vertical ellipses denote self- 
assembled chains of molecules, (b) A typical shape of a dipo- 
lar chain. Shaded ellipses denote single molecule wavefunc- 
tions in individual layers (see Eq. JTJ) with radius Rj. 



for dipoles located on top of each other in different layers. 
The system collapse in the direction of attractive inter- 
action is prevented by the strong optical potential. In 
the remainder of this paper we focus on electric dipolar- 
interactions between polar molecules but our conclusions 
should also apply to magnetic dipolar atoms [l[. 

Our main observation is that beyond a certain criti- 
cal value of the dipole moments, attractive interaction 
between molecules in different layers can bind them into 
chains. Longer chains appear first and have a stronger 
binding energy (see Fig. 2). At finite temperature there 
is a competition between entropy that favors shorter 
chains and energy which is minimized for longer chains. 
This leads to a non-monotonic dependence of the dis- 
tribution function on the length of the chains (see Fig. 
[3]). This behavior of the dipolar chains liquid (DCL) 
strongly resembles rheological electro- and magnetoflu- 
ids [9j]. In both cases the dipolar interaction leads to 
self-organization of elementary objects (nanoscale grains 
in rheological fluids and molecules in our case) into chains 
of varying length. In the rheological fluids several finite 
temperature phase transitions have been predicted the- 
oretically or numerically [To| . The chaining phenomena 
can also be directly imaged in a 2D thin film 11]. The 
new feature of DCL in polar molecule systems is the pos- 
sibility of Bose-Einstein condensation of chains at low 
temperatures. We find that the longest chains condense 
first and that in systems with many layers the conden- 
sation temperature is strongly suppressed compared to 
the transition temperature of individual molecules in one 
layer (insert to Fig. 3). Such suppression can be under- 
stood by observing that the presence of chains of differ- 
ent lengths reduces the number of dipoles available for 
the longest chains and hence limits the number of par- 
ticles available for the condensation. Finally we discuss 
implications of our results for experiments. 

We start by investigating the structure of a single chain 
in a multi-layered system in Fig. 1. For a strong ID 
optical lattice there is no spatial overlap of wavefunctions 
of molecules localized in different layers, hence M dipoles 
in successive layers can be described by the Hamiltonian: 
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Here Vj(r) - D 2 (|r| 2 - 2fd 2 ) / (|r| 2 + j 2 d 2 f' 2 is the 
dipole-dipole interaction with D being the electric dipo- 
lar moment in c.g.s. units, w is the in-plane trapping 
frequency and m is the mass of a single molecule. When 
L is the total number of layers in the ID optical lattice, 
we can have M = 2, 3, • • • , L as possible chain lengths. 
For simplicity we assume that all layers are equivalent 
and contain the same number of molecules. We define 
new variables: = Pj — P/M and Sj = Tj — R, where 

P = YljLi Pj an d R = jj YljLi r j are the total momen- 
tum and the center-of-mass(CM) position of the chain. 
As a result, Hm can be separated into sum of the CM 
part (H CM — j^. + |Afmw 2 R 2 ) and and the part that 
describes the relative motion inside the chain: 
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The commutation relation between the new coordinates 
are: [R, P] = ih, [sj, qj] = ihSji (l — M _1 ) , and all other 
commutators are zero. As we demonstrate below, im- 
portant properties of a DCL are dominated by the long 
chains (M ^> 1). So from now on we will treat Sj, as 
canonically conjugate variables for chains of all lengths. 
To calculate the binding energy from H ro i, we take the 
following variational wavefunction 
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where the wavefunction width, Rj , is determined by the 
competition of attractive part of dipole interaction, which 
favors localized wavefunctions, and the kinetic energy, 
which opposes localization. Near the chain center inter- 
action has contributions from both sides of the chain, 
while near the chain tips interaction comes from one side 
only and is weaker. Therefore we expect a smaller size 
of the wavefunction near the center compared to the tips 
of the chain (see Fig. D^b)). For simplicity we assume 
Rj = R (l + £|j - (M + 1)/2| 2 ), where R and £ are 
two variational parameters to be obtained by minimiz- 
ing (ipL({sj})\H Te \\%(jL({Sj})) . Taking more sophisticated 
trial wavefunctions will not have a qualitative effect on 
the properties of DCL that we discuss in this paper. 

In Fig. [5] we compare binding energies (per molecule) 
of the shortest (M = 2) and longest (M = 81) chains in 
a system with L — 81 layers. The control parameter is 
the dimensionless strength of the dipole interaction, Uq = 
mD 2 /h 2 d, which can be tuned either by increasing the 
electric field (which increases the value of dipole moments 
of molecules), by changing the interlayer distance d, or 
by increasing the effective mass of molecules using an 
additional in-plane optical lattice. A non-trivial result 
of our analysis is that bound states of many molecules 
(long chains) appear first and have a stronger binding 
energy. In the system under consideration, the longest 
chain (M = L = 81), can exist when Uq exceeds a critical 
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FIG. 2: Binding energy (note the minus sign) per molecule 
as a function of dipole strength, Uo, for chain length M = 2 
and M — 81. Results are obtained using the variational wave- 
function in Eq. Q, with u = O.OLEo. E = h 2 /md 2 = ^E R 
with Er being the recoil energy of the ID optical lattice. In- 
set: Optimal values of the variational parameters for a single 
chain, Rq (solid curve) and £ (dashed curve), as a function of 
Uo for M = 81. Note that the value of £ has been multiplied 
by 10 3 to fit the scale of the vertical axes. 



value Uq L ^ 0.4. A bound state of two molecules appears 
only forVo > C/ * 2 - 1.6 12]. So for U < C/q.l no 
chains are formed, for Uq > Uq 2 chains of all lengths 
are stable, and for intermediate dipole strengths (Uq l < 
Uq < Uq 2 ) there is a critical value of the chain length 
Mo such that chains of length Mo or larger are stable but 
shorter chains are not (except for the unbound dipoles, 
M = 1). In the inset of Fig. [51 we show the optimal 
values of the variational parameters, Rq and £, for the 
longest chain (M = 81) obtained as a function of Uq > 
Uq L . Our results show that when Uq is small (weakly 
bound chains), the minimum width of the chain can be 
several times larger than the interlayer distance, d. In 
this limit the curvature of the shape, £, is large and the 
chain has a shape of a horn in both ends (Fig. []Jb)). 
For large dipole strengths we find tightly bound chains 
in which both Rq and £ are very small. 

We now proceed to discuss the thermodynamics of 
DCLs. We start by making several simplifying assump- 
tions. Specifically, we neglect repulsive interactions be- 
tween chains and do not include bending of individ- 
ual chains. The first assumption is justified rigorously 
in the limit of small number of molecules within each 
layer. For N molecules in one layer, the characteris- 
tic strength of the in-plane dipolar energy is E- m t ~ 
ND 2 /a^ sc , where a osc = (h/muj) 1 ^ 2 is the in-plane oscil- 
lator length. We can neglect E- mt compared to hw when 
N 2 uj < (h 5 jD A m i ) 1 l 2 . We further note that most con- 
clusions involving thermodynamic properties of DCL will 
remain qualitatively similar even when interactions be- 
tween chains become important so long as they are not 
strong enough to lock molecules within the planes into 
a Wigner crystal phase @. Neglecting bending of in- 
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FIG. 3: Number of chains in a multilayer system (L = 81) 
as a function of the chain length, M. Other parameters are 
iVtot = 10 5 , Uo = 2, and huj = O.OLEo. Condensate T c is 
0.11-Eo for such parameters. Solid, dashed, and dotted lines 
are for temperature T = 0.1, 0.3, and 0.5 Eo respectively, and 
the upward arrow indicates the condensation of the longest 
chain (M — L — 81) Inset: Condensate T c as a function of 
dipole strength, Uo, for multilayer system L = 2 and L = 81. 
The number of dipoles per layer for these two systems are 
set the same so that they have the same T c in noninteracting 
regime. Solid lines are calculated within the noninteracting 
chain approximation (discussed in the text). Dotted lines are 
schematic since reliable theory in the vicinity of Uq l should 
include interactions between chains. 



dividual chains can be justified when temperatures are 
smaller than the characteristic binding energy and the 
chain lengths are not too large. This assumption may 
not necessarily apply near the condensation temperature 
for modest dipolar moments (see Fig. 3 and discussion 
below) but should be valid for large values of Uo [13J . 

When the interactions between chains are neglected, 
the total energy of a DCL is a sum of energies of individ- 
ual chains. Neglecting the bending modes, we can write 
the single chain energy as a sum of the center of mass 
and binding energies: em,\ — huj(i x + iy + 1) — E^m- 
Here i = (i x ,i y ) are quantum numbers for the CM co- 
ordinate in a 2D harmonic trapping potential and -E^m 
is the binding energy of a chain of length M. To un- 
derstand thermodynamics of the DCL we need to sum 
over all chain configurations with the constraint that the 
total number of dipolar molecules is fixed. We enforce 
the constraint using the Lagrange multiplier method and 
obtain he partition function, Z, 
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counts possible CM positions of chains in the z direc- 
tion. In the above partition function the chemical po- 
tential term Mfj, is proportional to the chain length M. 
This follows from the fact that the numbers of chains are 



not fixed but can change in dis/as-sociation type pro- 
cesses even in equilibrium. The only physical chemical 
potential fi is that of individual molecules and controls 
the distribution of chain lengths in thermal equilibrium. 
This is similar to well known systems of multi-ionized 
plasmas [l4[ and classical rheological fluids To find 
fj, in ([!]) we relate it to the total number of molecules 

iVtot = k B TdlnZ/d f i = Em=iMo MN **> where Nm 
is the number of chains of length M . At temperatures 
above Bose-Einstein condensation we have 
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Here zm = exp [f3(Eb,M + Mfi — Hu))] is the fugacity 
of chains of length M, and Lii(z) = X)fc°=i T? ^ s 
the polylogarithm function. In deriving expression ([3]) 
we replaced summation over discreet levels of the in- 
plane trap by the integration and used the identity, 

Jo°° Jo°° [ z ~ leX+y ~ l] dxdy = Li2{z). The prefactor, 
(L — M+l), counts the number of possible CM positions 
in the z direction. Bose-Einstein condensation means the 
appearance of a macroscopic occupation of some micro- 
scopic state. In our case this happens when the argu- 
ment inside the log of Eq. ^ goes to zero for some 
particular value of the quantum numbers M and i. From 
the dependence of the binding energy Eb,M on the chain 
length M, one can see that the condensation occurs for 
the longest chains, M = L, in the state i x = i y = 0. 
To find the transition temperature we need to solve the 
equation /i(T c ) = L~ l (hw — Ej, j,) = fi c . For tempera- 
tures below T c the chemical potential stays fixed at fj, c 
and the number of chains in the condensate is given by 

the equation N c = L^ 1 N tot - Y^m1\,m MN m (ii c ) . 

In Fig. [3] we show the calculated number of chains, 
Nm, as a function of their length for different tempera- 
tures. We choose Uo = 2.0 which is large enough to allow 
formation of the chains of all length (see Fig. [J). At high 
temperatures(dotted line), however, the entropy is more 
important, so short chains and unbound molecules dom- 
inate. One can show that in this limit the distribution 
function obeys the Saha's equation that was originally 
derived for multi-ionized plasma [l4| 
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As the temperature is lowered to be of the same order 
as the binding energy, the competition between entropy 
and energy leads to a nonmonotonic distribution of Nm 
(dashed line), which is a special feature of the DCL ap- 
proaching the the regime of quantum degeneracy. When 
the temperature is reduced further (solid line), quantum 
statistics becomes important and we find a strong en- 
hancement of the population of longer chains and even- 
tually Bose-Einstein condensation of the longest chains. 

In the inset of Fig. [31 we show the condensate T c 
as a function of the dipole strength for a multilayer 
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(L = 81) system. Solid lines in the Uq — > regime are 
obtained treating molecules as noninteracting particles. 
For large Uo regime we include self-assembly of molecules 
into chains but neglect interactions between chains. In 
the intermediate regime of Uq (around the critical dipole 
strength, Uq l ), we expect the approximation of nonin- 
teracting chains to break down due to a larger transverse 
size of the chains (see the inset of Fig. [5]). We use dot- 
ted lines to schematically describe the expected results 
of T c in this regime. We anticipate a sharp kink in the 
transition temperature in analogy to mixtures of bosonic 
atoms near Feshbach resonance regime [l5j], due to the 
formation of real space bound states (chains) at Uq l . We 
point out that T c is strongly suppressed in a multilayer 
system for Uo ~ Uq l because the number of the longest 
chains is very small. This is a special feature of dipolar 
chain liquid, since in this regime binding energies of the 
longest chain are small so that molecules are spread out 
among chains of various lengths. 

We now discuss experimental issues related to the real- 
ization and detection of the dipolar chain liquid proposed 
in this paper. For typical polar molecules with electric 
dipole moment D ~ 1 Debye, m ~ 100 atomic mass 
unit, and trapped in the optical lattice of Nd:YAG laser 
(d = A/2 = 0.53 ^irn), we find the dimensionless dipole 
strength, Uq = mD 2 /H 2 d, is about 2.8, and exceeds the 
critical dipole strength needed for the formation of dipo- 
lar chains (see Fig. [2]). Uq can be further controlled by 
an external DC electric field (Enc)- E.g. by starting 
from a BEC of the non-polarized molecules at Edc = 
one can drive a system into a DCL state by adiabatically 
increasing electric field. The resulting dipolar chain liq- 
uid can be observed by several different methods. For 
example, the frequencies of collective modes that are in- 
homogeneous along the stack will change dramatically 
for Uq > Uq L , due to the large change in the stiffness of 



the dipolar cloud associated with the emergence of long 
chains. Likewise, RF spectroscopy can be used to probe 
the energy gap resulting from the chain formation. Suc- 
cessive absorption peaks shall be observable in the exci- 
tation spectrum, yielding a signature of chaining effect in 
the dipolar cloud. Finally, a long-range density correla- 
tion between dipoles within the chains can be measured 
via noise correlation techniques after the cloud is released 
from the in-plane trap. 

We now discuss some novel avenues opened by this 
work. Firstly, starting with a fixed number of layers and 
approaching thermodynamic limit by increasing trans- 
verse system size, we expect that there will be a phase 
transition between a simple superfluid and a DCL at 
Uq = Uq L . Secondly, by increasing the dipole strength 
even further, a first order transition from superfluid of 
chains to a Wigner crystal phase should occur. In the 
Wigner crystal phase chains should form a triangular lat- 
tice with a fixed distance between them |6j . Finally, when 
considering bending and/or vibration of long chains at fi- 
nite temperatures, we expect the DCL to be similar to a 
classical fcrroliquid with nontrivial thermodynamic prop- 
erties. We expect that each of these avenues correspond 
to an intriguing direction for future theoretical and ex- 
perimental research. 

In summary, we have predicted the chaining effect of 
polar molecules confined in a ID optical lattice. The 
ground state and thermodynamic properties of the sys- 
tem can be qualitatively changed by the formation of 
chains. These intriguing phenomena should be observ- 
able in the experiments with ultra-cold polar molecules. 
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